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Abstract 

We prove a sufficient condition for synchronization for coupled one-dimensional maps and estimate 
the size of the window of parameters where synchronization takes place. It is shown that coupled sys¬ 
tems on graphs with positive eigenvalues (EVs) of the normalized graph Laplacian concentrated around 
1 are more amenable for synchronization. In the light of this condition, we review spectral properties 
of Cayley, quasirandom, power-law graphs, and expanders and relate them to synchronization of the 
corresponding networks. The analysis of synchronization on these graphs is illustrated with numerical 
experiments. The results of this paper highlight the advantages of random connectivity for synchroniza¬ 
tion of coupled chaotic dynamical systems. 
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1 Introduction 


Among problems concerning dynamics of large networks, synchronization holds a special place ll^ [33l . 
With applications ranging from the genesis of epilepsy |[T51 to stability of power grids ETl . understanding 
principles underlying synchronization in real world networks is of the utmost importance. From the theoret¬ 
ical standpoint, the mathematical analysis of synchronization provides valuable insights into the role of the 
network topology in shaping collective dynamics. 

The mechanism of synchronization in coupled dynamical systems in general depends on the type of 
dynamics generated by individual subsystems. For many diffusively coupled systems, the contribution of 
network organization to synchronization can be effectively described by the smallest positive eigenvalue 
of the graph Laplacian, also called the algebraic connectivity of the network ifTTll . This has been shown 
for coupled scalar differential equations (so-called, consensus protocols) ll30l 1^ . limit cycle oscillators 
1231, excitable systems forced by noise l25l . and certain slow-fast systems (261. Interestingly, stochastic 
stability of the synchronous state can be estimated in terms of the total effective resistance of the graph 
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12411371. Synchronization of the coupled chaotic systems features a new effect. It turns out that the parameter 
domain for synchronization depends on both the smallest and the largest positive eigenvalues of the graph 
Laplacian. This has been observed for systems of coupled differential equation models in ll^ |28l and for 
coupled map lattices ifT^ . Furthermore, the analysis in ifl^ shows that coupled systems on graphs, whose 
graph Laplacians have all positive EVs localized, are optimal for synchronization in the sense that they 
synchronize in a broader range of parameters compared to the systems on the graphs whose eigenvalues are 
spread out. The conclusions in lfT6l are based on formal linear stability analysis of synchronous solutions of 
coupled systems of chaotic maps. 

The goal of the present paper is twofold. First, we rigorously prove a sufficient condition for synchro¬ 
nization for systems of coupled maps. Our condition is slightly weaker than that proposed in HU, but it 
admits a simple proof and captures the mechanism of chaotic synchronization. Specifically, if shows fhaf 
diffusive coupling counferacfs infrinsic insfabilify of chaofic sysfem fo produce sfable spatially coherenf 
solufions. The synchronizing effecl of fhe coupling is more pronounced on graphs wifh localized posifive 
eigenvalues. To fhis end, in fhe second parf of fhis paper, we discuss specfral properties of cerfain symmefric 
and random graphs in fhe lighf of fhe sufficienf condifion for synchronizafion. Specifically, we review fhe 
relevanf fads abouf fhe eigenvalues of Cayley graphs I® and confrasf fheir properties fo fhose of quasir¬ 
andom and power-law graphs l|71|4l and expanders |[T3l . These resulfs highlighf fhe advanfages of random 
connecfivify for chaofic synchronizafion. 


2 Stability analysis 


Lef r = {V, E) be an undirecfed graph on n nodes. The sef of nodes is denofed by E = [n] := {1, 2,..., n}. 
The edge sef E confains (unordered) pairs of adjacenf nodes from V. Unless sfafed ofherwise, all graphs in 
fhis paper are assumed fo be simple meaning fhaf fhey do nol confain loops (ii ^ E Vi G V) and multiple 
edges. Denofe fhe neighborhood of f G U by 

N(i) = {jeV: ij G E}. 

The cardinalify of N{i) is called fhe degree of node f G U and denofed by di = |A^(i)|. If di = dVi £ [n] 
fhen r is called a d-regular graph. 


Af each node of T we place a dynamical system 

Xk+i = f{xk), ( 2 . 1 ) 

where / is a continuous funcfion from fhe unif inferval I := [0,1] fo ifself. Local dynamical systems af 
adjacenf nodes inferacf wifh each ofher via diffusive coupling. Thus, we have fhe following coupled sysfem 

xtli = f{Xk^) + j N>^ = 0,1,2,..., (2.2) 

* jeN(i) 

where e > 0 confrols fhe coupling sfrengfh. 


If r is a laffice, (2.2i is called a coupled map laffice. In fhis paper, T can be an arbifrary undirecfed 
connecfed graph. For convenience, we rewrife (|2.2[) in fhe vecfor form 


x^_i_i — iTff(x/j), Kg — Iji eZ/, 


(2.3) 
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where X = ..., and f(xfc) = (/(x^^^),/(x*^^^),...,/(x^"’^)). L stands for the normalized 

graph Laplacian of F: 

(2.4) 


L = h- D-^A 


where /„ is the nxn identity matrix, D = diag(di ,d 2 ,---, dn) is the degree matrix, and A is the adjacency 
matrix of F : 

1, ij E E, 

0, otherwise. 


= 


Thus, 


1, i = j, 

J-1 


{L)ij = { -d- , ij E E, 


0, otherwise. 

In general, L is not a symmetric matrix. However, it is similar to a symmetric matrix 

L = = In- D-^l‘^AD-^1’^. 

Thus, the EVs of L are real and nonnegative by the Gershgorin’s Theorem ifTil : 


(2.5) 


0 < Ai < A 2 < • • • < A„. (2.6) 

Furthermore, the row sums of L are equal to 0. Thus, Ai = 0 and 

V> = span{l„}, where := (1,1,..., 1 )t E M", (2.7) 

is the corresponding eigensubspace. If F is connected then 0 is a simple EV of L ifTTI . We will need the 
following properties of pseudosimilar transformations (cf. ll24ll F 

Lemma 2.1. Let L be the nxn Laplacian matrix of a connected graph F, whose EVs are listed in ( |2.6| ). 
Suppose S E is such that ker(S) = E and let S~^ denote the Moore-Penrose pseudoinverse of S 

mt. 

Then L = SLS^ is a unique solution of the matrix equation 

LS = SL. (2.8) 

The EVs of L counting multiplicities are 

A 2 , As,..., An- (2.9) 

The eigenspaces ofL corresponding to \i, i = 2,3,..., n, are mapped isomorphically to the corresponding 
eigenspaces of L by S. 

Proof The statement of the lemma follows from Eemma 2.4 of Il24ll . 


In the remainder of this paper, we use 



/ 
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5 = 


0 

-1 

1 .. 

. 0 

0 



( 2 . 10 ) 
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Note that S has full row rank and ker(5) = V. For a given x G one can use jS'xl to estimate the 
distance from x to V. Specifically, for the projection of x to V, Px>, we have 

|P 2 ?x| = |5+5x| < ||5+|||5x|. (2.11) 

Here and below, || • || stands for the operator norm fl4]l . 


V is an invariant subspace of (2.3 1 . Trajectories from this subspace correspond to the space homoge¬ 
neous solutions of (2.31. Thus, we refer to 2? as a synchronous subspace. We are interested in finding 
conditions on e, which guaranfee (asympfofic) sfabilify of V. In fhe analysis below, we follow fhe approach 
developed for sfudying synchronization in systems wifh continuous time 


Theorem 2.2. Let f : I ^ I be a twice continuously differentiable function and T = (H, E) be a connected 
graph on n nodes. Suppose 

F = max \ f{x)\ > 1, (2.12) 

and the EVs of the graph Laplacian ofT satisfy 

+ (2.13) 

'^min 

where Xmin cmd Xmax denote the smallest and largest positive EVs of the normalized graph Laplacian L, 
respectively. 


Let Xk, k = 0,1,2,... denote a trajectory o/(|2.2|) with 


ee{Xf\l-F-^),X-\l + F-^)). 


(2.14) 


Then there exists (5 > 0 such that iP^x^l 0 as k ^ oo provided \PtiXq\ < 5. 


For the proof this theorem, we will need the following auxiliary lemma, which we state first. 
Lemma 2.3. Let ^ sequence of nonnegative numbers such that 

k-l 

Uk < P^Uk-i A: = 1, 2, 3,..., 

j=0 


for some 0 < /r < 1 and A > 0. Then there exist positive numbers qo <^nd q such that 

Uk < qpf", k = 0,1,2,..., 


(2.15) 


(2.16) 


provided 0 < uq < qo. 


Proof. (Theorem |2.2| ) Let x^, A: = 0,1, 2,..., denote a trajectory of (2.31. By multiplying both sides of 

(2.17) 


(2.31 by S', we have 


yfc = SKff{xk) = K,Si{xk), K, = I^-i - eL, 
where = S'x^. Decompose 

Xfc = Ckln + Vfc, Ck = n~^llxk, Vfc = S'+S'xfc G V^, 


(2.18) 
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where we are using the fact that S~^S provides the projection on ker(S')-'“ = V-^. Using Taylor’s formula, 
we have 

S'f(xfc) = S {/(Cfc)l + /'(Cfc)vfc + gA:(vfc)} = /'(Cfc)yfc + gfc(yfc), (2.19) 

where y^. = S'xfc = 5vfc, gfc(yfc) = SgkiS+Yk), 

|gfc(y)l < C'l|y|^ (2.20) 


and Cl is a positive constant independent of A: G N U {0}. By plugging (2.191 in (2.17 1 , we have 

yfc+i = fkyk + gfc, 

where fk := f{ck) and gk := gfc(yfc). 

By iterating (|2.21[), we obtain 


( 2 . 21 ) 


Yk = 


liifA) 

j=0 


k-1 


yo + Yl 

j=0 


k-1 


n 

i=j+i 




( 2 . 22 ) 


where the product over an empty index set is assumed to be equal to the identity matrix. 

the EVs of 

L lie between A 2 and A^. Thus, fi := F||I„_i — eL|| < 1, whenever 

eG(A2-i(l-F-'),A-ni + i"-')). 


Condition (2.131 implies that (A2^(1 — T’ ^),A„^(1 + F ^))is nonempty. By Lemma 


2.1 


Thus, 


|/ill|i^e|| < /i < 1, j G Nu {0}. 


(2.23) 


Using (2.201 and (2.23), from (2.22) we obtain 


fc-i 


\yk\ < fi^lyol + Cl k eN. 

j=0 


Using Lemma [23| from (2.24) we obtain 


|yfc| < C 2 ^i^, A: G NU {0}, 


provided yo is sufficiently small. 

□ 


Proof. (Lemma [23] ) 

We use induction to prove ( 2.16| ). Let ?? > 0 be fixed. Denote 

/(O = ^ exp{-A^//i(l - n)] 


and set ryo := /(r?) < rj. Clearly, (2.16) is true for A = 0 and 0 < uo < rj. 


(2.24) 


(2.25) 
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Suppose (2.161 holds for 0 < fc < p for some p G N. Below we show that this implies that it also holds 
for k = p.To this end, let 

, k — 0 , 1 , 2 ,.... 


Using (2.281, we have 


p-i 

Zp ^ Uq + A. ^ ^ fi Up—l—jZp—l—j 
j=0 

p-1 


< 


uo + Arj'^ pP ^ ‘^Zp-i-j 


j=0 

p 


= uo + Ar]'^p^ 
i=i 

where we used the induction hypothesis in the second inequality. 

Using the Gronwall’s inequality (see Lemma [2^ below), from (2.261 we obtain 

p 

Zp < no exp{^p ^ < noexp{A?7/p(l - p)} = 


< r], 


i=i 


(2.26) 


(2.27) 


where we used 0 < no < p < /(p) in the last inequality. Thus, 

Up < r]pP. 

The statement of the lemma follows by induction. 

□ 

The proof of following discrete counterpart of the Gronwall’s inequality can be found in IITtII . 
Lemma 2.4. Let {zk}^Q and {pfcj^i be two nonnegative sequences such that 


Zk<B + fe G [p], 

i=i 


(2.28) 


for B > 0 and p G N. Then 


Zk < B exp Pj 
i=i 


> , fc G [p]. 


3 Connectivity 


In this section, we examine the role of the network topology for chaotic synchronization in the light of Theo¬ 
rem!; 


[ 2.2 The sufficient condition for synchronization (2.131 favors graphs whose nonzero EVs are localized. 
Thus, we review several families of graphs r„ = {V (L^), i7(r„)), n G N, which possess this property even 
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Figure 1: a) A Cayley graph based on a ball, b) a Paley graph, and c) an ER graph. 


when |l/(r„)| grows without bound as re —)• oo. These families include Erdds-Renyi (ER) random, quasir¬ 
andom, and power-law graphs ciiiiia, and expanders ifTSll . We also review spectral properties of Cayley 
graphs. Because Theorem [2]^ provides only a sufficient condition for synchronization, we tested how well 
(2.14 1 captures the domain of synchronization numerically. We found a good agreement between numerical 
results and the analytical estimates of the synchronization domain. Numerics also confirm fhaf localization 
of fhe graph EVs is critical for synchronizafion of large nefworks. 


3.1 Cayley graphs 

We sfarf our review of nefwork connecfivify wifh Cayley graphs. These are highly symmefric graphs defined 
on groups. 

Definition 3.1. Suppose G is a finite additive group and S C G is a symmetric subset, i.e., S = —S. Let 
r = {V, E) be a graph defined as follows V = G and for any a,b€G,ab€Eifb — a£S. T is called a 
Cayley graph and is denoted Cay(r, S). 


The following Cayley graphs on fhe (addifive) cyclic group 
connecfivify in (2.2 1 . 


Z/reZ, re G N, model nearesf-neighbor 


Example 3.2. Let n > 3,r < [re/2j and B{r) = {±1, ±2,..., ±r}. T = Cay(Z„, i?(r)) is called a 
Cayley graph based on the ball B{r). 


Cayley graph T = Cay(Z„, 5(1)) corresponds fo fhe sfandard nearesf-neighbor inferacfions in fhe 
coupled system (2.2). By varying r in Cay(Z„, B{r)), one can control the range of coupling (see Pig|^). 
The analysis of the coupled system (2.2 1 on T = Cay(Zn,, B{r)) below shows that for fixed r and re 1 


fhey do nof meef synchronizafion condifion ( 2.13| ). This will be confrasfed wifh fhe analysis of nefworks on 
random and quasirandom graphs, which will be discussed in §3.2| 

The EVs of Cayley graphs on Z„ can be compufed explicifly in ferms of fhe one-dimensional represen- 
fafions of Z„ : 

= (ej;(0),e2;(l),.. .,Gx{n - 1))^, x G Z„, (3.1) 
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where 


exiy) = exp I I , ?/ E Z„. 


For the EVs of Cayley graphs on a cyclic group, the following lemma is well-known ifT^ . 
Lemma 3.3. The EVs of the normalized graph Laplacian ofV = Cay(Z,i, S) are given by 

Aa; = 1 - ^ea;(s) = 1 - |S|"^ ^COS , X eZn. 

y&S yes V / 

The corresponding eigenvectors are e^, x E Z„. 


(3.2) 


(3.3) 


Proof. For any x,y € 1in,'we have 


{Ae^){y) = ^x{y + s) = X] 


ses 


^ses 


Thus, 


AGx - hx^X) hx - ^ ^ Sj; (s) ■ 

seS 


(3.4) 


Since characters e^;, x E Z„, are mutually orthogonal, (3.4) gives the full spectrum of A. The statement of 


the lemma follows from (3.4i and L = |S'|/„ — A. 

□ 

Corollary 3.4. The spectrum o/T = Cay(Z„, S) contains a simple zero EV Aq = 0 corresponding to the 
constant eigenvector !„ = (1,1,..., 1)^. 


With the explicit formula for the EVs of the Cayley graph Cay(Z„, S) in hand (cf. (3.31), we study syn¬ 
chronization in the coupled system ( |2.2| ) on symmetric graphs. We start with the nearest-neighbor coupling 
corresponding to E = Cay(Z„, 


(3.5) 


Example 3.5. By applying ( |3.3| ) toT = Cay(Zn, 5(1)), we obtain the following formula for the nonzero 
EVs of F: 

(2'kx\ 

Ax = 1 — cos - , X E n — 1 . 

\ n J 


By setting x = 1 and x = in (3.5), we find the smallest and the largest positive EVs ofF, respectively: 

= 1 - COS ( — 1 = -h 0{rr‘^), (3.6) 

(3.7) 


A. 


n J 

2, if n is even, 

1 — cos(7r(l — n~^)) = 2 + 0{n~‘^), ifn is odd, 


Thus, 


A 


max ft —2\ 

= +0{n ). 


(3.8) 


We conclude that the nearest-neighbor family of graphs does not satisfy synchronization condition ( |2.i3| )/Qr 
large n. 













Figure 2: The range of the variance of a typical trajectory of (2.2 1 on Cay(Z„, i?(l)) as a function of the 
coupling strength e: a) n = 5, b) n = 6, and c) n = 7. 


Following ifT^ . in our numerical experiments we use the following measure of coherence. For {x^, k G 


N}, a trajectory of (2.2 1 , consider 


{a\x))k = -Xfc)2, xfc = = 0,1,2,.... 


(3.9) 


i=i 


i=i 


We will refer to <t^(x) as the variance of x. Clearly, synchronization means that |((T^(x))fc| stays small 
for A: >> 1. Numerical results for (2.2 1 on several nearest-neighbor graphs are shown in Fig. The plots 
in Fig. show the range of values of (after removing transients) for a typical trajectory for the range 
of coupling strength e G (0,1] for the nearest-neighbor graphs of different size. The windows in e where 
the range of is close to 0 indicate synchronization. Plot a) for n = 5 shows a window of synchrony 
for an interval in e, which shrinks for n = 6 (Fig. |^b); and already for n = 7 (Fig. |^c) no windows of 
synchronization were found. 

We now turn to F = Cay(Z„, B{r)) for r G N and n > 2r. This is the case of the r-nearest-neighbor 
coupling. 


Example 3.6. For T = Cay(Z„, B{r)), using ( |3.3| ), we have 

Xx = 1 


1 


f 2'Kxk\ 


\k=—r 


= 1 — 



— 1 I , x G [n — 1]. 


(3.10) 


Using the Taylor’s formula, after straightforward calculations from {3.10) we estimate 

7r^(2r -h l)(r -h 1) 


0 ^ Xrnin ^ -^1 — 


‘in? 


+ 0{n 


-4\ 


(3.11) 


To estimate Xmax use x* = [ 2 ^+ 1 ]- For fixed r > 1 and n ^ 1, using Taylor’s formula and elementary 
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Figure 3: The range of the variance of a typical trajectory of (2.21 on Cay(Z„, B{r)) as a function of the 
coupling strength e: a) n = 100, r = 10, b) n = 200, r = 20, and c) n = 400, r = 40, d) n = 100, r = 25, 
e) n = 200, r = 50, and f) n = 400, r = 100, g) n = 100, r = 49, h) n = 200, r = 98, and i) n = 400, 


r = 196. 
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inequalities, we estimate 


sin 


TT \ / TTX \ ( TT 

- < Sin - < sin - 

2r+ly \ n J \,2r + l 


—7r(2r + 1) 


+ 0{n < sin 


n 


7r(2r + l)x* 


n 


< sin 


Thus, 

and 


Xr 


> 


3n^ 


Xmin 7r^27’(r + 1) 


+ 0{n 




+ 0(n ^), 


TT 


2r + 1 


Xmax ^ Xx* > 1 + (2r) + 0{n ) 


+ 0(n 


(3.12) 


Estimate (3.12 1 shows that for any fixed r G N the coupled system (|2.2|) may satisfy (2.131 for small 


n and fails to do so starting with some critical network size n*{r) (see Fig.^a-c). For larger values of r, 


satisfies (2.13 1 for larger n. In fhis respecf, fhe presenf case is nof differenf from fhe nearesf-neighbor 
coupled nefworks, which we discussed in Example |3.5| The sifuafion changes if we consider nonlocal 
coupling, i.e., r = [pn\ wifh fhe range of coupling p G (0,0.5). Then clearly for nof too small p, one can 
have ( |2.13| ) even torn S> 1. 


3.2 The complete, random, and quasirandom graphs 


The complefe graph on n nodes = {V, E) is defined hy V = [n] and E = {ij : I < i < j < n}. 

Lemma 3.7. The normalized graph Laplacian of has a simple zero EV and the EV 1 + (n — 1)“^ o/ 
multiplicity n — 1. 


Proof. Since Kn is a connected (n — l)-regular graph, fhe zero EV of L{Kn) is simple, 
corresponding eigenvector (cf. (2.7 1 ). Nexf, note fhaf any nonzero vector v = {vi,V 2 , ■ ■ ■, Vn) 
In is an eigenvector of L{Kn). Indeed, suppose Vi Q for some i G [n]. Then 


and In is fhe 
orfhogonal fo 


{L{Kn)v)i = Vi- 



1 


n — 1 


where fhe orfhogonalify condifion = 0 used. 

□ 


Therefore, fhe synchronizafion condifion (2.131 always holds for T = Kn, n > 2. Furthermore, for 
sufficienfly large n, fhe synchronizafion interval ( |2.14| ) is nonempfy, and is pracfically independenf of n. 
The windows of synchronizafion shown in Fig.[^a-c do nof change as fhe size of fhe graph increases from 
n = 10 in Fig. fo n = 100 in Fig. [^. Below we consider ( |2.2| ) on several graphs fhaf approximate Kn 
and, like Kn, exhibif good synchronizafion properties. 

We begin wifh fhe classical Erdos-Renyi (ER) model of random graphs ||6l. In fhis model, fhe edge 
befween every pair ij, 1 < z < j < n, is inserfed wifh probabilify p. The decision for a given pair is made 
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Figure 4: The range of the variance of typical trajectories of (2.2i on (a-c), G{n, 1/2) (d-f), and Paley 
(g-i) graphs plotted versus the coupling strength parameter e. The graph sizes are a) n = 10, b) n = 50, c) 
n = 100, d) n = 10, e) n = 50, f) n = 100, g) n = 13, h) n = 53, and i) n = 101. 
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Figure 5: The range of the variance of typical trajectories of (2.2 1 on the power-law graphs on 200 nodes. 
The degree distribution of the graphs used in all three plots is ~ k~^. The average degree is a) d = 20 and 
h)d = 30. 


independently from the decisions on other pairs. The resultant graph on n nodes is denoted by G{n,p) (see 
Fig.[^). The nonzero EVs of G{n,p) satisfy 

irmx|l — AjI < (3.13) 

almost surely (cf. |[9ll3T]ll. Therefore, from the synchronization viewpoint ER random graphs are as good as 
complete graphs. They also have considerably fewer edges than Kn, because the expected degree of edges 
in G{n,p) is pn{n — l)/2. 

The EV concentration results in fQ] apply to some other random graph models including power-law 
graphs (see also ma). Specifically, let G{w) be a random graph on n nodes with a specified expected degree 
sequence w = {wi,W 2 , ■ ■ ■, Wn) f8|. Chung-Eu random graphs on a vertex set V = [n] can be constructed 
by inserting an edge between i and j from V with probability pij = WiWj {Ylk=i Here, we assume 

maxjgy wf < Ylk=i '^k and allow for loops. Eurther, for a given power law exponent /3, maximum degree 
m, and average degree d, one can generate a power-law graph by setting (cf. ifTOll l 

Wi = c{i + 1 < f < rr, 


where 


c = 


/ 9-2 

/3-1 


dnV/3 


-1 


and fo = 


d{P - 2 ) 

m{f3 — 1) 


0-1 


Eor Chung-Eu power-law graphs. Theorem 9 in ifTOl then guarantees that all positive EVs are localized 


|Afc(L)-l| =0(1), k>2, 

with probability 1 — o(l) provided minjgy Wi S> In n. In Eig.j^we present numerical results for the coupled 
system ( |2.2| ) on Chung-Eu power-law graphs for two values of the expected degree, d. Eor both values of d, 
Eig.j^shows large windows of synchronization. 

Next, we turn to quasirandom graphs, which share many combinatorial properties with the ER graphs 
G{n,p) (cf. niHl). Below, we restrict to the case p = 1/2 and consider only d-regular graphs. There 
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are several equivalent properties that define a quasirandom graph T = {E,V) Q. For instance, for any 
S' C y, {\V\ = n), the number of edges connecting the nodes in S is 


-\S\^ + 0{n\ 


(3.14) 


exactly as for the ER graph G(n, 1/2). If F is a d-regular graph then (3.14i is equivalent to the following 
localization property of the nonzero EVs of the normailized graph Eaplacian (cf. Theorem 9.3.1, Q) 


|Ai - 1| = o(l), z = 2,3,... ,n. 


(3.15) 


Hence, for coupled systems on quasirandom graphs synchronization condition (2.131 holds for sufficiently 


large n and, as in the case of the complete graph, the interval of synchronization (2.14i to leading order is 
independent of re forn » 1 (see Eig. -f). To illustrate synchronization properties of ( |2.2| ) on quasirandom 
graphs, we use Paley graphs, which are defined as follows. 

Example 3.8. Let re = 1 (mod 4) be a prime and denote 

= Z„/{0} and Qn = {x^ (mod re) : x G Z^}. 

Qn is viewed as a set (not multiset, i.e., each element has multiplicity 1). Then Qn is a symmetric subset of 
Zf and \Qn\ = 2“^(re — 1) (cf. (|7S| Lemma 7.22]). Pn = Cay(Z„, Qn) is called a Paley graph l[18V . 


The nonzero EVs of L{Pn) can be computed explicitly from (3.3), using properties of the Gauss sum 
(cf mi): 

X (p\ = I ^ ^ ^ n 161 

^1 n) I xisnotaQR (mod re), 


where x G Z^. 


The numerical resulfs for fhe coupled sysfem (|2.2|) on Paley graphs show good agreemenf wifh fhe 


analyfical esfimafes. The synchronizafion windows compufed numerically for (2.21 on Paley graphs are as 
large as fhose for fhe complefe graphs of fhe same size (see Pig.|^g-i). Overall, fhe variance plofs for (2.2i 
on fhe complefe, ER, and Paley graphs in Eig.j^have a greaf degree of similarify. 


3.3 Expanders 


The complefe, random, and quasirandom graphs endow fhe coupled sysfem (2.2 1 wifh good synchronizafion 
properfies buf af a price of having 0(n‘^) edges in a graph of size re. In fhis subsecfion, we will review cerfain 
families of graphs wifh much fewer edges, which nonefheless approximafe complefe graphs and promofe 


synchronizafion in coupled sysfems like (2.2). 


Eef = {V(Tn), E{Tn)), re G N, be a family of d-regular connected graphs. Suppose V (F^) = [re] 
and denote fhe EVs of L(Tn) by 

0 = Ai < A 2 < • • • < Xn- 
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By Gershgorin’s Theorem Ifl4l . < 2d. Therefore, the most likely way for the synchronization condition 
(2.131 to break down for large n is through the second EV A 2 (r„) approaching 0. This leads us to consider 
families of graphs, for which the second EV remains bounded away from zero 


A2(rn) > Cl > 0. 


(3.17) 


Such families of graphs are called expanders 1351. 

Definition 3.9. /[3?117^/ A family of d-regulargraphs E^ = (E(rn), i7(rn) (|V(r,i)| —)• oo <35 re —>■ ooj is 
called an expander family if the Cheeger constant 


h(Tn) = min^--— ■ ■ j ——-——• 


> C 2 > 0 Vre, 


(3.18) 


where the minimum is taken over all proper subsets ofV{Tn), F and dF stands for the set of edges 
connecting the nodes in F with the nodes in the complement of F, V (r„)/i^. 


The Cheeger constant h(r) quantifies the connectivity of T. The uniform bound on /i(r„) in (3.181 
guarantees that r„ remain well-connected even as |V(r„)| grows without a bound. The combinatorial 


condition (3.18) is equivalent to the spectral one (3.17) (cf. mmi which we arrived at in the context 
of synchronization. Therefore, synchronizability in large networks is directly related to connectivity: the 
better connected large networks are, the more likely they to satisfy the synchronization condition ( |2.13| ). 
The connectivity is better in graphs with larger second EV A 2 ifTTl . There is an upper bound on A 2 (r„) 
(cf. HI): 

limsupA2(r„) <d — 2\/d — 1. (3.19) 

n^oo 

A family of graphs {T^} is called Ramanujan if 


A2(rn) > d — 2y/d — 1. 


Ramanujan graphs exhibit the optimal asymptotics of A 2 (r„) and, in this sense, possess the best possible 
connectivity. By inspecting several common families of d-regular graphs, such as the nearest-neighbor 
graphs, one can see that typically the second EV tends to 0 rapidly as re —)• 00 (see, e.g., ( |3.11| )). In fact, any 
family of Cayley graphs of bounded degree on an abelian group can not be an expander family ifT^ §4.3]. 
Nontheless, a typical family of random d-regular graphs (d > 3) is very close to be Ramanujan with high 
probability as shown in ifT^ . The deterministic constructions of expander families are also available, albeit 
they rely on sophisticated algebraic techniques (cf., EHEIlEl). 


To illustrate synchronization in coupled systems on expanders, we will use the following family of 
random bipartite d-regular graphs. 

Example 3.10. Let B 2 m,d, m > 2 be a bipartite graph on 2m vertices. The edges are generated using the 
following algorithm: 


1. Let p : [m] —)• [rre] be a random permutation. In our numerical experiments, we used MATLAB 
function randperm to generate random permutations. For i G [m], add edge (i, rre -|- p{i)). 

2. Repeat step i. d — 1 times. 
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Figure 6: The range of the variance of typical trajectories of (2.2i on the random bipartite graphs (cf. Ex- 
ample [3T0| . The graphs used in these simulations have 400 nodes. The degrees are a) 10, b) 15, and c) 20. 
For comparison, in plots d-f, we included numerical results for the coupled system on Cayley graphs with 
regular connections Cay(Z,i, B{r)) for d) n = 100, r = 10, e) n = 200, r = 20, and t) n = 400, r = 40. 



Figure 7: The range of the variance of typical trajectories of (2.2i on the random Cayley graph 
Cay(Z„, ±Rk): a) n = 100, A: = 10, b) n = 200, k = 20, and c) n = 400, k = 40. 
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The numerical results in Fig. show that coupled system on i? 4 oo,i 5 already has a notable window 
of synchronization (see Fig. [^), which becomes even bigger for i? 4 oo ,20 (see Fig. 15). Note that in the 
present case we do not have localization results for the EVs of L, as in the analysis of the coupled system 
( |2.2| ) on random and quasirandom graphs (cf. ( 3.13| ), ( |3.15| )), and we can only rely on the lower bound for 
A 2 ( |3.17 1 . Therefore, for chaotic dynamical systems on expanders, in general, we can not guarantee that 
synchronization condition (2.131 holds for all n » 1. Nonetheless, as can be seen from the numerical 
results in Fig. [^a-c, one can achieve synchronization on large expanders of relatively small degree. The 
numerical results for the bipartite random graphs may be compared to the results for the Cayley graphs on 
B{r) in Fig. |^d-f, which show that the coupled system (2.2) on smaller graphs of the same degree but 
with regular connections remains rather far from synchrony. This is another illustration of the advantages of 
random connectivity for synchronization. 


In conclusion, we consider synchronization of the coupled system (2.2 1 on random Cayley graphs 
Cay(Z„, where Rj- is a set of k elements chosen from the uniform distribution on ■ Random 

Cayley graphs are expanders with high probability forn » 1 and k = 0(log n) |[20ll . The numerical results 
for the dynamical systems on random Cayley graphs in Fig. [7] show that these systems exhibit much better 
synchronization properties than the systems on regular Cayley graphs of the same size (see Fig.[^d-f). 


4 Discussion 


Motivated by the linear stability analysis in iTT^ . in this paper, we rigorously derived a sufficient condition 
for synchronization in systems of diffusively coupled maps, whose dynamics can be chaotic. In contrast to 
the approach in ifT^ . we did not seek to show asymptotic stability of spatially homogeneous solutions of 

dTH ) 

Xfc = Xkln, k = 0,1,2,, (4.1) 


where {x^} is a trajectory of the local system ( 2.1[ ). Instead, we proved the asymptotic stability of the 
invariant subspace of such solutions V (cf. (2.7 1 ), in the sense that every trajectory of (2.2) that started 
sufficiently close from V approaches V in forward time. Note that the asymptotic stability of V does not 
imply the asymptotic stability of individual trajectories ( |4.1| ). In particular, our approach circumvents the 
analysis of perturbations in the tangential to V direction, which would be necessary for the justification of the 
linearization about ( |4.1[ ). The latter problem is quite technical already in the case when local systems have 
stable dynamics If23ll . On the other hand, it is not clear that for coupled chaotic systems one has sufficient 
control over perturbations in the tangential direction, because it corresponds to a degenerate subspace of the 
coupling matrix and the local dynamics are sensitive to perturbations. 


The sufficient condition for synchronization ( 2.13| ) identifies the complete graph as an optimal network 
organization for chaotic synchronization. Therefore, in the second part of this paper, we considered networks 
on graphs that approximate the complete graph. These include ER random graphs, as well as quasirandom 
and Chung-Eu power-law graphs. We show that from synchronization point of view these graphs are almost 
as good as the complete graph provided that the network is large enough. Eurthermore, expanders, which 
constitute a broader class of graphs, tend to have good synchronization properties. As many of these graphs 
are constructed using random algorithms, our results highlight the advantages of random connectivity for 
synchronization. 
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